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A theoretical analysis of the statistical distributions of the reflected intensities from random media 
is presented. We use random matrix theory to analytically deduce the probability densities in the 
localization regime. Numerical calculations of the coupling to backward modes in surface corrugated 
waveguides are also put forward for comparison. Interestingly, the speckle distributions are found 
to be independent of the transport regime. Despite the scattering being highly non-isotropic, the 
predicted probability densities reproduce accurately the numerical results. 



PACS numbers: 42.25.Bs, 41.20.Jb, 5.45. +b, 84.40.Az 



The scattering of both classical and quantum waves 
from random media has long been a subject of interest 
PJ. Interesting multiple scattering effects like the en- 
hanced backscattering (EB) and intensity correla- 
tions in transmitted and reflected waves have exten- 
sively been analyzed. Recently much attention has been 
focused on this subject in connection with the probability 
distributions of the different transmittances in disordered 
waveguides [0dTlJ ■ Rayleigh and Gaussian statistics with 
anomalous tails in the diffusive regime, observed exper- 
imentally [0], were analytically deduced within the ran- 
dom matrix theory (RMT) framework A strong 
dependence on the length L of the disordered region was 
expected from numerical simulations of waveguides with 
volume disorder |1J. Recent results ||, also based on 
RMT, have shown how the Rayleigh and Gaussian dis- 
tributions evolve into the same lognormal distribution 
as L increases beyond the localization length £. The 
same behavior has been found in simulations of surface 
corrugated waveguides However, in contrast with 
the transmission case, the probability distributions of re- 
flected intensities have not been discussed in detail pre- 
viously. It is the purpose of this work to put forward a 
comprehensive picture of the statistics of reflected waves 
from a disordered waveguide based on analytical as well 
as on numerical calculations. 

In order to discuss some general properties of the back- 
ward scattering, random matrix theory provides a use- 
ful macroscopic approach since it does not depend on 
the specific details of the scattering processes. We will 
make use of RMT results to analytically deduce the prob- 



ability density of the reflected speckle pattern in a N- 
mode waveguide in the localized regime. As expected 
from the central limit theorem, in the limit of large 
N, the speckle distribution follows a negative exponen- 
tial law (Rayleigh statistics). As a consequence of en- 
hanced coherent backscattering effects (i.e. time-reversal 
invariance S), as N decreases, this Rayleigh law evolves 
into two different distributions depending on the back- 
ward scattering directions. The numerical calculations 
of the speckle statistics for a corrugated waveguide show 
a remarkable agreement with these analytical results all 
the way from the diffusive through the strong localized 
regime. Although the analytical results are formally de- 
duced in localization, this agreement indicates that, in 
contrast with transmission, the reflected speckle distribu- 
tions are independent of the transport regime. 

Our theoretical approach is based on the properties 
of the scattering matrix S of a waveguide with a cen- 
tral region containig the random medium (whose disor- 
der can be indistinctly either in volume or in the walls) 
(see Fig.l). Outside this region, the transverse confine- 
ment in the waveguide defines propagating modes or 
channels |Q. The wave transport can be described in 
terms of the scattering matrix S: 
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where r, t, t, f, are square matrices whose dimension 
is given by the total number of channels N. The ma- 
trix elements r,j and U m axe the reflected and transmit- 
ted amplitudes into the channels j and to, respectively, 



1 



when there is a unit flux in the channel i incident from 
the left; r nm and t n j have the same meaning but the 
incoming channel is incident from the right. In terms 
of these elements, the transport coefficients Ru — \rij\ 2 



and TL 



= U 



measure the reflected and transmit- 



ted speckle patterns respectively, being Ri — J2j Rij an d 
Tjj = ^2 rn Ti m the total reflection and transmission for 
the incoming mode i. 

The study of the probability densities of the reflected 
intensities P(Rij) is much simpler in the strong localiza- 
tion regime where the transmission amplitudes vanish. 
The matrix S is then reduced to the two square blocks r 
and f. Since we are interested in a statistical approach, 
given our limited knowledge of the microscopic coupling 
between channels, the natural choice of the statistical en- 
semble is that which maximize the information entropy 
subject to the known constrains [ flux conservation (S 
unitary) and time-reversal invariance (S symmetric) ]. 
In the RMT context, this leads to the circular orthogo- 
nal ensemble (COE, f3 — 1), i.e. the phases of the ele- 
ments of the reflection blocks are as random as possible 
while their moduli are constrained by the conditions of 
unitarity and symmetry. In the absence of time reversal 
invariance, the matrix S loses its symmetry and the ap- 
propriate statistical ensemble is then the circular unitary 
ensemble (CUE, (3 = 2). 

Following general formulas for the distribution of ma- 
trix elements in the circular ensembles derived by Pereyra 
and Mello jl5j , we can compute the distribution of the 
reflected speckle pattern P(Rij). Due to the underlying 
isotropy condition p 12 the RMT results are mode in- 
dependent except for the backscattering into the same 
channel: 

P(R*,P = 1) = (^-y - l) (1 - R it ) 1/{R " } - 2 , (2) 
and 
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where 2F1 is the hypergeometric function and the mean 
values {Rij) are given by 



1 + SjjS/n 
N + Sf3i 



(4) 



Although recent numerical calculations |ll|,|16|,^7[ show 
that the different channels are not equivalent, as we will 
see, the actual statistical distributions can still be de- 
scribed by the general RMT results. The key point will 
be to consider the distributions as functions of the mean 
values rather than in terms of the number of channels. 

In Figs. 2 and 3 we have plotted Eqs. (Q) and (||) for a 
set of values of (Ra) and (Rij), respectively (thick solid 



and broken lines), together with numerical results for a 
surface corrugated waveguide (see the details below) . To 
this end, (Rij) always involves i 7^ j. In the limit of N 
large , (both (Ra) and (Rij) < 1), both distributions (||) 
and (||) evolve into a Rayleigh law. Coherent backscat- 
tering effects manifest themselves in the distributions for 
small N which dramatically depend on whether i — j, or 
i ^ j. This is illustrated in the analysis of the case N = 2 
((Ra) = 2/3 and (Rij) = 1/3) where the probabilities for 
i = j and for i ^ j are quite different (see inset in Fig. 
3): 



and 



P(Riul3=l) 



P(Rij,(3 = l) = 



2VT^R~, 



(5) 



(6) 



When time reversal symmetry is removed, there is no 
longer difference between P(Rij) and P(Ru), and hence 
there is a unique intensity distribution. Interestingly, in 
terms of the averages the statistical law coincides with 
that obtained for Ra in the case of {3 = 1, i.e. 

P(R i j,0 = 2)=P(R i i,p = 2) = P(R u ,0=l) . (7) 

This means that the statistical law that governs the 
backscattered intensities Ra does not depend on the time 
reversal symmetry conditions, as long as they are ex- 
pressed in terms of the mean value. 

It is worth noticing, that for a chaotic cavity |l8||, Eqs. 
U and H are exact results tor the reflected intensities for 
any N. In this case, (Rij) = (1 + S lj S /31 ){2N + tfpi) -1 . 
Again, in terms of the mean values, the statistics of re- 
flected waves from a chaotic cavity and from a waveguide 
in the localization regime are indistinguishable. 

Let us now discuss the numerical results for the reflec- 
tion coefficients in a surface disordered waveguide (see 
Fig. 1). The corrugated part of the waveguide, of total 
length L and perfectly reflecting walls, is composed of n 
slices of length I. The width of each slice has random 
values uniformly distributed between Wq — 5 and Wq + 5 
about a mean value Wo- Following previous works, we 
shall take Wo/ 8 = 7 and 1/8 = 3/2. The main trans- 
port properties do not depend on the particular choice 
of these parameters, however ]lT| , [l9| . Transmission and 
reflection coefficients are exactly calculated by solving 
the 2D wave equation by mode matching at each slice, 
together with a generalized scattering matrix technique 
]l9|-pT[. In order to obtain the intensity distributions as 
well as the mean values for the reflection coefficients, we 
have performed the calculations over an ensemble of one 
thousand configurations. 

The anisotropy of the scattering process, predicted in 
Refs. [ pTp^ , p^ ], is specially relevant when the EB factor 
r)i, defined as 
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r )i = (N- 1) 



(Ri) — (Ra 



(8) 



is addressed. This EB factor is very sensitive to the 
changes in the mean values, and so it is specially affected 
by changes in the incident wavelength. In Fig. 4b we 
show the behavior of 77 as a function of the wavelength, 
for a length of the corrugated part such that the system 
is always in the localization regime. The behavior of the 
"external" modes (defined as those propagating modes 
with either the smallest or the largest transversal mo- 
mentum) is quite different to that of the "central" ones. 
The lowest mode (i.e. that with the smallest transversal 
momentum) has a marked oscillatory value of 77, show- 
ing peaks at the onset of each new propagating mode 
in the waveguide. For the rest of the modes, except for 
small oscillations, 77 evolve from a high value just at the 
onset of propagation, to the expected factor of 2. This 
oscillating behavior is closely related to that observed in 
the localization length £ shown in Fig. 4a, (see also Ref. 

The numerical probability densities have been calcu- 
lated in the range of wavelengths shown in Fig. 4, 
for different length values including both the diffusive 
(£/N < L < £) and localized (L > £) regimes. The 
histograms in Fig. 2 and 3 (thin solid lines) show the 
intensity distributions numerically obtained for different 
averages of either (Ra) or (Rij). Distributions corre- 
sponding to different number of propagating modes and 
different transport regime, but having the same average 
of the reflection coefficient, are indistinguishable. Despite 
being the scattering highly non-isotropic, it is remarkable 
that both numerical and analytical results coincide with 
each other, without any adjustable parameter. This in- 
dicates that, while the average values depend on the spe- 
cific system under consideration, the fluctuations, which 
are characteristics of disordered systems, are fully rep- 
resented by the probability density as a funtion of the 
mean values, and have general properties that are well 
described by RMT. 

In summary, we have derived the general statistical 
properties of the reflected intensities from a random 
medium in the localization regime, based on RMT. Com- 
parisons have been made by numerically analizing the 
coupling to backward modes in surface corrugated waveg- 
uides. The enhanced backscattering factor can be much 
larger than that predicted by RMT, and oscillates ver- 
sus the wavelength following the oscillating behavior ob- 
served for the localization length. We show, however, 
that in spite of the anisotropy of the interactions, the 
numerical statistical distributions accurately follow the 
RMT predictions, independently of the transport regime. 
This indicates that the statistical distributions do not de- 
pend on the details of the scattering processes. They 
should then be observable in different systems under 
rather general disorder conditions. 
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FIG. 1. Schematic view of the system under consideration. 




FIG. 2. Distributions of Ru for different values of (Ru). 
Thin solid lines correspond to the numerical calculations for 
(Ru) ~ 0.15,0.51, respectively, whereas thick solid lines are 
the analytical predictions. Thick broken lines are the analyt- 
ical results for (Ru) ~ 0.25,0.35 and 0.45, respectively. 
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FIG. 3. Distributions of Rij(i 7^ j) for different values 
of (Rij). Thick solid lines are the analytical predictions for 
the same values of the averages, as those indicated for the 
numerical histograms (thin solid lines) {{Rij} ~ 0.05,0.2, re- 
spectively). Thick broken lines are the analytical results for 
(Rij) ~ 0.10 and 0.15, respectively. Inset: Probability density 
for N = 2: P(R ii ) ({Ru) = 2/3), and P(R lJ ) ((R tj ) = 1/3). 
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FIG. 4. (a) Localization length £ versus Wo/A. (b) 
Enhanced backscattering factor r)i versus Wo/A. Vertical 
long-dashed lines are drawn at the onset of a new mode. 
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